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ABSTRACT 

Two series of low-cycle-fatigue (LCF) test data for two different groups of aircraft gas turbine 
engine compressor disk geometries were reanalyzed and compared by using Weibull statistics. Both 
groups of disks were manufactured from titanium (Ti-6A1-4V) alloy. A probabilistic computer code 
called Probable Cause was used to predict disk life and reliability. A material-life factor A was 
determined for the titanium (Ti-6A1-4V ) alloy by using fatigue disk data and was applied to predict 
disk life as a function of speed. A comparison was made with results from the currently used life 
prediction method, which is based on crack growth rate. Applying an endurance limit to the computer 
code did not significantly affect the predicted lives under engine operating conditions. Failure 
location predictions correlated with experimentally observed failure locations in the LCF tests. A 
reasonable correlation was obtained between the disk lives predicted by using the Probable Cause 
code and those predicted by using a modified crack growth method. Both methods slightly 
overpredicted life for one disk group and significantly underpredicted it tor the other. 


NOMENCLATURE 

A material-life factor 
c stress-life exponent 
e Weibull slope 

F probability of failure, fraction or percent 
h exponent 
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constant, deg C (deg F) 


K T 

L life, number of stress cycles or hr 

Lq | 0.1 -percent life, or life at which 99.9 percent of a population survives, number of stress 

cycles or hr 

L| 1 -percent life, or life at which 99 percent of a population survives, number of stress cycles 
or hr 

L| 0 10-percent life, or life at which 90 percent of a population survives, number of stress 
cycles or hr 

L 5() 50-percent life, or life at which 50 percent of a population survives, number of stress 
cycles or hr 

Lp characteristic life, or life at which 63.2 percent of a population fails, number of stress 
cycles or hr 

m Weibull slope 

N life, number of stress cycles 

n number of components, elemental volumes, or bolt holes 
S probability of survival, fraction or percent 
T steady-state temperature, °C (°F) 

AT thermal gradient, deg C (deg F) 

V stressed volume, (in.-^) 

X load, time, or stress 

X|, X-y fractional time at load and/or speed 
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Z depth to maximum critical shear stress, m (in.) 
a stress or strength, N/m 2 (psi) 
o u location parameter, N/m 2 (psi) 

x critical shear stress, N/m 2 (psi) 

x 45 maximum shear stress, N/m 2 (psi) 

Subscripts 

/ ith component or stressed volume 
n number of components or elemental volumes 
ref reference point, stress, volume, or life 
sys system or component probability of survival or life 
v related to stressed volume 

P designates characteristic life or stress 


INTRODUCTION 

According to Coy , 1 a simplistic approach to component life and reliability would be a completely 
deterministic one. That is, each component in the system would be subjected to the scrutiny of 
determining its calculated stress for an assumed load in service. The calculated stress would be 
compared with an allowable stress for an assumed material strength. The deterministic approach 
assumes that if the material strength is greater than the applied stress, the component will not fail (i.e., 
the life will be infinite). (Stress is compared with strength by using the appropriate failure criteria, 
e.g., maximum shear stress, strain energy, distortion energy, or fatigue.) The deterministic method 
assumes that full and certain knowledge exists for the service conditions and the material strength. 

For most aerospace applications designing for infinite life with a zero probability of component 
failure would require structural size and weight that would not be compatible with aircraft and 
spacecraft flight. As a result it becomes prudent to design for finite life with an acceptable risk or 
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probability of failure. Hence, a more reasonable assessment of component life and reliability for 
aerospace and, perhaps, nonaerospace applications is based on a probabilistic approach to design and 
analysis. 

A. Paimgren 2 - 7 in 1924 suggested a probabilistic approach to predicting the lives of machine 
components and, more specifically, rolling-element bearings. On the basis of his test results he 
suggested that an acceptable life be defined as that time at which 10 percent of a population of 
bearings will have failed or 90 percent will have survived. He also noted that there was an apparent 
size effect on life. That is, larger bearings with the same equivalent load as smaller bearings had 
shorter lives than the smaller bearings. 

In 1939 W. Weibull 4 ' 6 published two papers that describe a statistical approach to determining 
the strength of solids. Weibull postulated that the dispersion in material strength fora homogeneous 
group of test specimens can be expressed according to the following relation: 



( 1 ) 


where X = a and Xp = Gp (see Appendix A). 

Equation (1) relates the probability of survival 5 and the fracture (or rupture) strength G. When 
In In ( 1/5) is used as the ordinate and In g as the abscissa, fracture (and fatigue) data are assumed to 
plot as a straight line. The slope (tangent) of this line is referred to as the “Weibull slope” or “Weibull 
modulus” usually designated by the letter e or m. The plot itself is referred to as a “Weibull plot.” 

By using a Weibull plot it becomes possible to estimate a cumulative distribution of an infinite 
population from an extremely small sample size. The Weibull slope is indicative of the dispersion 
of the data and its density (statistical) distribution. Weibull slopes of 1 , 2, and 3.57 are indicative of 
exponential, Raleigh, and normal (Gaussian) distributions, respectively. 7 

The scatter in the data is inversely proportional to the Weibull slope. That is, the lower the value 
of the Weibull slope, the larger the scatter in the data and vice versa. The Weibull slope is also liable 
to statistical variation depending on the sample size (data base) making up the distribution. 8 The 
smaller the sample size, the greater the statistical variation in the slope. 

Weibull 4,5 related the material strength with the volume of the material subjected to stress. If we 
imagine the solid to be divided in an arbitrary manner into n volume elements, the probability of 
survival for the entire solid can be obtained by multiplying the individual survivabilities together as 
follows: 


5 = S r 5 2 ...5„ 


( 2 ) 


where the probability of failure F is 


F=\-S 


(3) 
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Weibull 4 ’ 5 further related the probability of survival S, the material strength o, and the stressed 
volume V, according to the following relation: 


"4 = J v /awv 


(4) 


where 


f(X) = G e 


(5) 


For a given probability of survival S 


o ~ 



( 6 ) 


From Eq. (6) for the same probability of survival the components with the larger stressed volume 
will have the lower strength (or shorter life). This principle was applied to successfully normalize 
rolling-element fatigue life by Carter 9 and Zaretsky et al. 1() This equation was also applied by 
Grisaffe 1 1 in a Weibull analysis of the bond strength of plasma- sprayed alumina coatings on stainless 
steel. Grisaffe showed that the calculated mean bond shear strength decreased with increasing test 
area in accordance with Eq. (6). 

In 1947, G. Lundberg and A. Palmgren 12 applied Weibull analysis to the prediction of rolling- 
element bearing fatigue life. The Lundberg-Palmgren theory expressed /(X) in Eq. (4) as 


/(*) = 


x c N e 

Z h 


(7) 


where x is the critical shear stress, N is the number of stress cycles to failure, and Z is the depth to 
the maximum critical shear stress in a concentrated (Flertzian) contact. From Eqs. (4) and (7) 


N ~ 



Me 


( 8 ) 


From Lundberg and Palmgren* 2 (see Appendix B) the lives of individual stressed volumes at a 
given probability of survival are summarized as follows: 



( 9 ) 
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Equations ( 1 ) to (4) and (7) to (9) are the basis for the international standards for rating rolling- 
element bearing life. 15 

In Lundberg and Palmgren 12 the parameter c/e is the stress-life exponent. This implies that the 
inverse relation of life with stress is a function of the life scatter or data dispersion. A search of the 
literature for a wide variety of materials would suggest that the stress-life exponent is independent 
of Weibull slope e. Also, from observation, and contrary to the Lundberg-Palmgren assumption, 
fatigue life appears to be independent of the depth to the maximum critical shear stress in a body. 
Hence, Zaretsky 16 has modified the Lundberg-Palmgren theory where 

f(X) = X ce N e (10) 


From Eqs. (4) and (10) 


N ~ 



(ID 


Zaretsky 1 7 as well as loannides and Harris 1 8 proposed a generalized Weibuil-based methodology 
for structural life prediction that uses a discrete-stressed- volume approach. August and Zaretsky’ 9 
extended this methodology by developing a technique for predicting component life and survivabil- 
ity that is based on finite-element stress analysis. Holland et al. 20 applied this method to predicting 
the fracture strength and life of a metal-matrix composite ring and used coupon data to determine the 
Weibull parameters for the material. Melis and Ogonek 21 implemented this life prediction method- 
ology through a computer code called Probable Cause. 

The computer code Probable Cause has not been benchmarked and compared with a statistical 
data base for aircraft engine components. Such a data base does exist. Mahorter et al. 22 reported low- 
cycle-fatigue (LCF) testing conducted by the Naval Air Development Center (NADC) and the Naval 
Air Propulsion Center (NAPC) and the resultant life data for two different groups of gas turbine 
engine compressor disks. The disks were each manufactured from a titanium (Ti-6A1-4V) alloy. 
They compared the results with analytical predictions by a crack growth method. Their conclusion 
was that their analytical predictions were not conservative. 

In view of the aforementioned it became the objectives of the work reported herein: ( 1 ) to exercise 
the computer code Probable Cause to predict the life and reliability of the two groups of gas turbine 
engine compressor disks tested by the Navy; (2) to compare and benchmark the predicted and 
experimental results; (3) to determine a material-life factor for use in the computer code for a titanium 
(Ti-6A1-4V) alloy; and (4) to determine what effect, if any, the assumption of a material endurance 
limit may have on predicted disk life. 


COMPUTER CODE PROBABLE CAUSE 

The computer code Probable Cause is written in PATRAN command language (PCL) and is 
based on the computational method of August and Zaretsky. 19 The component life and survivability 
are predicted by finite-element analysis. By establishing a reference life L re( , a unit or gage volume 
V ret -, and a reference stress x related to the reference life and gage volume, Eq. (11) can be rewritten 
substituting L for N as follows: 
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L — Lfct' 


( 12 ) 


' tret 7 

TKcfl 

X 

. v _ 


By replacing X in Eq. ( 1 ) with L, the probability of survival S and the life L can be related to *ref and 
L re | as follows: 


JL rc( /L) e 

d 


(13) 


By using Eqs. (12) and ( 13) and critical shear stress results from a finite-element analysis, L and 
S values for each element are computed. Hence, the probability of survival for the entire analysis 
model can be obtained by using Eq. (2) to multiply the individual survivabilities. 

These equations provide relative or normalized values for L and S in relation to reference values 
chosen from the selected element. Generally, reference values of 1 .0 and 0.9 are assigned to the L ret 
and S re f variables, respectively, in the equations. These values imply a relative or normalized life of 
unity and a probability of survival of 90 percent for the reference element or volume V rcl . Maximum 
shear stress or Von Mises stress values are normally used because they are applicable to failure theory 
for ductile materials. A reference element or volume can be chosen at random. However, the authors 
have primarily used the element with the highest resultant stress as a reference point. The stress-life 
exponent c and the Weibull slope e are parameters specific to the material. August and Zaretsky 19 
have previously discussed the rationale for obtaining these values. 

The computer code was embedded in a PATRAN finite-element code. The program operates on 
analysis results within the PATRAN data base. It stores the life and survivability predictions for each 
stressed volume (mesh element) as additional result files for postprocessing by using the standard 
PATRAN tools. 2 1 The data output provides the system life and the probability of survival. In addition 
to the conventional stress contours for the component, contour plots are provided for probability of 
survival (or failure) and life. 


TURBINE DISKS 
Specimens and Procedure 

Mahorter et al. 22 studied two different groups of aircraft gas turbine engine compressor disks 
designated disks A and disks B. The disks were manufactured from a titanium (Ti-6A1-4V ) alloy. The 
material parameters for this alloy, taken from various published sources, are given in Table 1. 
Schematic drawings of these disks are shown in Fig. 1. Five unfailed disks of each type were 
randomly selected from those removed from engines at overhaul. The disk operating times are 
summarized in Table 2. These failure data are the only such data reported and available in the open 
literature. Using them allowed us to exercise the Probable Cause code and to compare the results with 
actual components that had been run in aircraft engine service. 

Data were not made available by Maborter et al. 22 as to the total number of unfailed disks of each 
type removed from service and their respective operating times, the number ot disks that were 
removed from service because of crack initiation and their respective operating times, and the 
material processing and metallurgical characteristics of the titanium alloy used for each disk. Hence, 
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we do not know the actual population size that each set of disks represents. Also not reported were 
the manufacturing process and engine operating temperature for each disk. These unreported 
variables can affect disk life and the accuracy of the statistical life results reported herein. 

In the aircraft engine these disks operated at speeds between 1 1 000 and 1 1 800 rpm. The disks’ 
in-flight operating temperatures were not reported. After removal at engine overhaul these disks were 
isothermally cyclic tested at a speed of 1 1 200 rpm in spin pits at the Naval Air Propulsion Center 
(NAPC ). 22 A spin pit test cycle is shown in Fig. 2. One cycle is equivalent to 1.45 flight hours for 
disks A and 1.54 flight hours for disks B. The equivalent stress cycles for engine operation and for 
the NAPC (spin pit) tests are given in Table 2. A failure is defined as the first occurring crack equal 
to or greater than 0.7938 mm (0.03 1 25 in.) no matter where in the disk it occurs . 22 These data were 
analyzed by using the method of Johnson 8 and two-parameter Weibull analysis because of the limited 
number of data points. Three-parameter Weibull analysis would estimate a running time below which 
no failure would be expected to occur ( 100 percent probability of survival). (See Appendix A.) The 
results of these analyses are shown in the Weibull plots of Fig. 3 and summarized in Table 2. 

For aircraft engine applications and comparisons a probability of survival of 99.9 percent is used 
(i.e., the operating time exceeded by 99.9 percent of the components or the time at which 0. 1 percent 
have failed). This time is designated as the L 0 ] life. According to Mahorter et al. 22 at the L 0 i life 
disk retirement policy recommends either the removal or reworking of the “used disk.” From Eq. ( 1 ), 
where X = L, the Lq j life is calculated from the data of Fig. 3 according to the following: 


A). I =^exp- 


6.9073 

e 


(14) 


Stress Analysis and Life Prediction 

A N ASTRAN linear static finite-element analysis (FEA) was used to predict the maximum shear 
stresses T 45 of the disks. The FEA model mesh of sections of disks A and B are shown in Fig. 4. 
Symmetry in the disks allowed for reduced model segments to be used for analysis. Hence, disk A 
was analyzed as a 60° segment and disk B as a 30° segment. 

Both disks were analyzed for stress and probability of failure for two load conditions at 1 1 200 
rpm. One condition related to disk operation in the engine under full compressor blade (rim) loading. 
The second condition related to the NAPC spin tests without compressor blade (rim) loading. These 
FEA maximum shear stress results for the disks are shown in Fig. 5. For these analyses the contact 
(Hertz) stresses between the bolt and the bolt hole and between the blade root and the blade slot were 
not considered. 

For both disks a common gage or reference volume was chosen with a reference stress. Although 
this volume and stress can be chosen arbitrarily, it is most efficient for calculation purposes to select 
a gage volume having a high stress. For purposes of these calculations a reference volume of 
2.264709x10“ 10 m 2 (1.3820105xl0“ 5 in. 2 ) having a x 45 stress of 509.2 MPa (73.846 ksi) was 
chosen. 

A sensitivity analysis performed by Zaretsky et al . 22 oh a generic rotating disk showed that the 
L l0 life was insensitive to changes in Weibull slope e greater than 2 and to stress-life exponents c 
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between 7 and 1 3. The results of these analyses are shown in Fig. 6. From these analyses and the data 
of Fig. 3 a Weibull slope of 2 was chosen for the current life analysis to represent the failure 
distribution of both disks. By using fatigue data from the Aerospace Materials Handbook 24 we 
calculated a stress-life exponent of 9.2 for the titanium (Ti-6A1-4V) alloy and assumed it to represent 
that of the disk. 

The results of these analyses showing the probability of failure and the failure location on each 
disk are shown in Fig. 7. This figure suggests that for both disks and at all load conditions the bolt 
holes have the highest probability of failure. All the failure locations reported by Mahorter et al. 22 
for both sets of disks were bolt hole cracks. According to Mahorter et al. 22 turbine engine LCF 
experience has shown that the bolt holes are the primary location of engine disk failure. 


RESULTS AND DISCUSSION 
Comparison of Lives 

By using the method of Zaretsky 16,17,19 and the computer code Probable Cause, 21 a life analysis 
was undertaken of two different groups of aircraft gas turbine engine compressor disks for which 
there existed fatigue data. 22 These disks were manufactured from titanium (Ti-6A1-4V) alloy. The 
normalized predicted lives are summarized in Table 3(a). The disks may have been operated under 
two load conditions. Life prediction for a disk under combined loading depends on the percentage 
of the total operating time at each load and the individual life attainable under that load. When using 
a linear damage (Palmgren-Langer-Miner) 2 ' 25 ' 26 rule, the predicted disk life under the combined load 
conditions is as follows: 


4ys L \ h. 

Referring to Table 2, the time at each load condition was varied for both disks A and B. However, 
it can be reasonably assumed, for purposes of calculation, that for disks A the time share was 18 
percent under engine load conditions and 82 percent under spin pit conditions and that for disks B 
the time share was 15 and 85 percent, respectively. The results of these combined load-cycle 
calculations are shown in Table 3. Comparing these results qualitatively would suggest that disks A 
would provide a life approximately 14 percent longer than that provided by disks B. 

According to Mahorter et al. 22 "the test buildup for the compressor disks was essentially a 
subsection of the compressor complete with spacers and blade.” It could therefore be reasonably 
assumed that the spin test loading approached or equaled the engine loading. As a result further 
comparisons were made for the loaded (engine) conditions. Under these conditions the lives ot disks A 
would be expected to be 35 percent longer than those of disks B. 

From Fig. 3 and Table 2 theL () , and L |() lives of disks B were longer than those of disks A, but 
the reverse was true for the L ^ life. Unfortunately, the number of data points for both disks A and 
B was not sufficient to provide statistical significance to the life differences, to the absolute life 
values, or to the failure distributions (Weibull slopes) shown in Fig. 3. 
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Effect of Endurance (Fatigue) Limit 


An endurance (fatigue) limit for a material suggests that where the stresses in the material are 
below this limit the life of that elemental volume will be infinite (i.e., the probability of survival is 
100 percent). If the entire structure has its stresses below this value, the structure would not be 
anticipated to fail from fatigue. From the literature 24 the estimated fatigue limit for the maximum 
shear stress x 45 for the titanium alloy used is 258.5 MPa (37.5 ksi). This value is derived from tensile 
loading fatigue tests. 

Referring back to the FEA stress analysis in Fig. 5, those elemental stressed volumes whose 
stresses are less than 258.5 MPa (37.5 ksi) were assigned a survival probability of 100 percent. The 
disk lives were recalculated accordingly. These life results are summarized in Table 3(b) for the two 
load conditions of each group of disks and were combined by using Eq. (15). For the unloaded 
conditions the fatigue limit can significantly affect life. However, for the loaded and combined 
conditions the assumption of a fatigue limit did not significantly alter the calculated lives summa- 
rized in Table 3(a). 


Reference Life and Material-Life Factor 

By using Eq. (12) the reference life L ref was normalized as one cycle for 
V ref =2.264709x10“ 10 m 3 ( 1.3820 105x1 0" 5 in. 3 ) and x ref = 509.2 MPa (73.846 ksi). In order to 
assign an absolute value to L re , to predict life, it is necessary to benchmark the normalized predicted 
lives in Table 3(a) with the experimental L j0 lives in Table 2. Since disks A resulted in the shortest 
experimental life, they were used to benchmark the analysis. 

For combined loading of disks A 


0.33147 L,. ef = 7900 cycles 
Lfgf = 23 833 cycles 

For the engine (loaded) condition of disks A 

0.05968 L p e f = 7900 cycles 

= 132 373 cycles 


(16a) 


(16b) 


A material-life factor A for the titanium (Ti-6A1-4V) alloy can be defined and determined from 
Eq. (12) where 

^ = A’ef ^ref ^ref (17a) 

For L ref from Eq. (16b), e = 2, c = 9.2, V ref = 2.264709x10“ 10 m 3 (1.3820105xl0 -3 in. 3 ), and 
x ref = 509.2 MPa (73.846 ksi) 
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A = 2.527 x 10 8<) cycles 


N 


m 


9.2 ' 

TK9 


(International units) 


(17b) 


A = 3.024 xlO 47 cycles- 
(Customary units) 

From Eq. ( 1 ) for a Weibull slope e of 2 the L () , life can be calculated as a function of the L |() life as 
follows: 


lb 


9.2 


tn. 


16.9 


^,=0.07447^0 (18) 

By using the reference lives from Eqs. ( 16a) and ( 16b) and the computer code Probable Cause, the 
absolute predicted L, 0 lives in number of cycles were calculated for disks A and B at four speeds. 
From the L, () live the L () , lives were calculated from Eq. (18) and converted into engine hours. 

These results are summarized in Table 4. The predicted t lives from Table 4(b) are plotted as 
a function of speed in Fig. 8 and compared with the experimental L () , lives from Table 2. At the lower 
speeds disks B are predicted to have longer life than disks A. At the higher speeds the reverse is 
predicted. This crossover appears to occur at or about 1 1 200 rpm. Had the Von Mises stress values 
been used for this analysis the resultant lives may have been altered from those predicted by using 
the maximum shear stress x 45 . 


Thermal and Residual Stress Effects 

The analysis reported was undertaken only with those material data and operating conditions 
available in the open literature. Because the disk operating temperatures were not reported, 22 we did 
not consider them. That is not to imply that temperature has no effect. The temperature and resulting 
temperature gradients in the disk have a two-fold effect. First, the fatigue life of a material generally 
decreases with temperature. Second, thermal gradients in the disk induce thermal stresses. 

Zaretsky et al. 2 ** present an expression for thermal effects (based on uniform, radial, linear 
temperature gradients) on disk relative life as follows: 


L ~ 


K t 

AT 


-,0.52 


(19) 


where AT is the total radial temperature difference and Kj is a constant equal to 0.56 deg C when 
ATis in degree Celsius or 1 deg F when ATis in degrees Fahrenheit and where ATis equal to or greater 
than Kj. Hence, if the temperature gradient in the disk is known, Eq. (19) can be used to adjust the 
life prediction up or down. 
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The bulk operating disk temperature will also affect life. In general, compressor disk operating 
temperatures can approach or exceed 649 °C (1200 °F). There are no available data for the disk 
material showing the effect of bulk temperature on fatigue life. 

Another variable not considered and for which information was not available was compressive 
residual stresses present in the disk material. These residual stresses could be present as the result of 
either heat treatment or shot peening. Accordingly, the compressive residual stresses can reduce the 
effective shear stresses and thus increase life. These data are also not available. 

The effects of temperature, thermal stresses, and residual stresses on the accuracy of the life 
predictions and the determination of the material-life factor A may or may not be significant to the 
work reported herein. Since for approximately 15 to 20 percent of the time the disks were run at 
elevated temperatures in the engine, the value of A as it applies to ambient (room) temperature 
conditions may be too low. However, at elevated conditions the value of A may be to high. 

The effect of residual stresses may be relatively the same at both room temperature and elevated 
temperatures. Thermal stresses were not present at room temperature in the spin pit tests. At ambient 
(room) temperature the residual stresses will increase life relative to the same material without 
residual stresses. 


Crack Growth Life Prediction Method 

The currently used analytical life prediction method is based on material crack propagation rates 
at the highest stress point in the structure. The time for an initial flaw size to reach a defined length 
or size determines the mean or average life of the structure. 

Mahorter et al. 22 did not report what this initial flaw size was assumed to be but defined the initial 
crack length to failure as being equal to or greater than 0.7938 mm (0.03125 in.). The logic in 
choosing this crack length is that it is the smallest that can be detected by current state-of-the-art 
inspection techniques. 

It is assumed that the density distribution about the log of the mean life prediction is normal (or 
Gaussian). Accordingly, a 3a deviation about the log of the mean life results in an L {) , (99.9 percent 
probability of survival) life prediction. 

For the compressor disk the bolt holes are the highly stressed area and are the location where 
Mahorter et al. 22 assumed crack initiation to occur. They detected no cracks in the blade slot area 
of disks A. Mahorter et al. 22 calculated disk life for only a single bolt hole location and did not 
consider multiple bolt hole locations. Their calculated L 0 f lives are summarized in Table 5. We 
modified the NAPC-predicted lives 22 by using Eq. (9) to account for the entire number of bolt holes 
in each disk, where Eq. (9) (see Appendix B) can be written as follows: 




( 20 ) 


and where L is the NAPC-predicted life for disks with a single bolt hole, n is the number of bolt holes, 
and e is the Weibull slope of 2. For disks A and B, n equal 12. The recalculated lives based on 
Eq. (20) for disks A and B are summarized in Table 5. Comparing the L 0 | life predictions by the 
modified current method with those by the Zaretsky Probable Cause method shows a reasonably 
close correlation. Both methods slightly overpredicted life for disks A and significantly underpredicted 
life for disks B. 
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SUMMARY OF RESULTS 


Two series of low-cycle-fatigue (LCF) test data for two different groups of aircraft gas turbine 
engine compressor disks were reanalyzed and compared by using Weibull statistics. Both groups of 
disks, designated disks A and disks B, were manufactured from a titanium (Ti-6A1-4V) alloy. Tests 
had been run by the U.S. Naval Air Propulsion Center (NAPC). A probabilistic computer code based 
on the methods of Zaretsky and called Probable Cause was used to predict disk life and reliability. 
Probable Cause, which is a method for estimating a component’s design survivability and life, 
incorporates finite-element analysis and probabilistic material-life properties. The analysis was 
normalized to one set of disk data to obtain a material-life factor A applicable to the titanium alloy. 
This material-life factor A was used with the computer code to predict disk life as a function of speed. 
A comparison was made with a modified current life prediction method, which is based on crack 
growth rate. The following results were obtained: 

1 . A reasonable correlation was obtained between the disk lives predicted by using the computer 
code Probable Cause and those predicted by using a modified crack growth life prediction method. 
Both methods sightly overpredicted the life for disks A and significantly underpredicted life tor 
disks B. 

2. A material-life factor A based on disk fatigue data was determined for the titanium (Ti-6 Al- 
4V) alloy and applied to predict disk life as a function of speed. The value of A for the titanium alloy 
was 2.527xl0 80 cycles • N 9 - 2 /m l6y (3.024xl0 47 cycles • lb 9 - 2 /in. 16 - 9 ). 

3. The failure probability analysis for both disk groups predicted with reasonable engineering 
certainty the failure locations at the bolt holes of each disk. These locations correlated to those 
observed experimentally in the disk LCF tests. 

4. Applying an endurance (fatigue) limit to the computer code Probable Cause did not 
significantly affect the predicted disk lives under engine operating conditions relative to when an 
endurance limit was not used in the analysis. 
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APPENDIX A— DERIVATION OF WEIBULL DISTRIBUTION FUNCTION 


According to Weibull 6 any distribution function can be written as 


F(X) = l-exp-[/(X)] (Al) 

where F(X) is the probability of an event (failure) occurring. Conversely, from the above the 
probability of an event not occurring (survival) can be written as 

l-F(X) = exp-[/(X)] (A2a) 

or 

1 -F = exp-[/(X)] (A2b) 

where F = F(X) and (1 - F) = 5, the probability of survival. 

If we have n independent components, each with a probability of the event (failure) not occurring 
being ( 1 - F), the probability of the event not occurring in the combined total of all components can 
be expressed from Eq. (A2b) as 


F n — exp — [n/(X)] (A3) 

Equation (A3) gives the appropriate mathematical expression for the principle of the weakest link 
in a chain or, more generally, for the size effect on failures in solids. As an example of the application 
of Eq. (A3), we assume a chain consisting of several links. Also, we assume that by testing we find 
the probability of failure F at any load X applied to a “single” link. If we want to find the probability 
of failure F„ of a chain consisting of n links, we must assume that if one link has failed the whole 
chain fails. In other words, if any single part of a component fails, the whole component has failed. 
Accordingly, the probability of nonfailure of the chain (1 - F ;( ), is equal to the probability of the 
simultaneous nonfailure of all the links. Thus, 

\-F n =(\-Ff (A4a) 


or 


S„ =S" (A4b) 

Or, where the probabilities of failure (or survival) of each link are not necessarily equal (i.e., 
S| * St ^ S 3 ...), Eq. (A4b) can be expressed as 
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S'! = iS’l • S2 • S3 • . . . 


(A4c) 


This is the same as Eq. (2) of the main text. 

From Eq. (A3) for a uniform distribution of stresses throughout a volume V 


F y - 1 - exp -[Y/'(c)] 
or 

S = 1 - F v = exp ~[Vf(a)] 
Equation ( A5b) can be expressed as follows: 


In In 



In /(o) + In V 


(A 5 a) 


(A5b) 


(A6) 


It follows that if In ln( 1/S) is plotted as an ordinate and In /(a) as an abscissa in a system of 
rectangular coordinates, a variation of volume V of the test specimen will imply only a parallel 
displacement but no deformation of the distribution function. Weibull 4 assumed the form 


/(C) 


a -a, 


til 


( A7) 


and Eq. (A6) becomes 


In In 


1 

S 


= m In ( o - G„ ) - m In CTp + In V 


(A8) 


If g„, which is the location parameter, is assumed to be zero and Vis normalized whereby In Vis zero, 
Eq. (A8) can be written as 


lnln 


' 1 ' 

= /?/ln 

a 


.5. 


_°p_ 


(A9) 


Equation (A9) is identical to Eq. ( 1 ) of the main text. 

The form of Eq. ( A9) where o M is assumed to be zero is referred to as “two-parameter Weibull.” 
Where a tl is not assumed to be zero, the form of the equation is referred to as “three-parameter 
Weibull.” The Weibull slope or modulus m, which is also designated by e, is the slope or tangent of 
the line when the equation is plotted on rectangular coordinates. 
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APPENDIX B— DERIVATION OF SYSTEM LIFE EQUATION 


G. Lundberg and A. Palmgren 12 in 1947, using the Weibull equation for rolling-element bearing 
life analysis, first derived the relationship between individual component lives and system life. The 
following derivation is based on but is not identical to the Lundberg-Palmgren 12 analysis. 

From Appendix A, Eq. (A9), the Weibull equation can be written as 



1 


N 

In In 


= m In 


S’ 

^sys 


! 

1 


(Bl) 


where N is the number of cycles to failure. 

Referring to the sketch of a Weibull plot in Fig. 9, the slope m can be defined as follows: 



1 


1 

Inin 

- Inin 


1 

C/5 

>> 

Co 5 


_^ref _ 


In N - In N k{ 


(B2a) 


or 


In 

1 

Skys. 


1 

.5 ref . 

In 


N 


-|/M 


From Eqs. (Bl) and (B2b) 


(B2b) 


and 


1 


1 

N ’ 

m 

N 


= 

In 

— 


1 

C/5 

1 


^ref _ 



v p _ 


(B3) 


S = exp- 


m 


JL 

Vp 


(B4) 


Referring to Fig. 10, for a given time or life N, each component or stressed volume in a system 
will have a different reliability S. From Eq. (A4c) for a series reliability system 


Ssys -^1 S 2 St, ... 


(B5) 
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Combining Eqs. (B4) and (B5) gives 


exp- 


m 


N_ 



N 

in 

N 

111 

n! 

exp- 

1 

+ exp- 

n 

sf - 
1 

+ exp- 

1 

1 


(B6a) 


or 


exp- 


N 

N a 


= exp- 


1 m r 


N 


N, 


PI 


N 


Nt 


P2 


N 


Nt 


P3 


in 


+ ... 


(B6b) 


It is assumed that the Weibull slope m is the same for all components. From Eq. (B6b) 


1- 1 

m 

- 

m 

- 

m 

- 

m 

N 

, 

N 

+ 

N 

+ 

N 

+ ... 

_N$_ 


3._ 




_^ P 3_ 



(B7a) 


Factoring out N from Eq. (B7a) gives 


N a 


Nt 


Pi 


m 


+ 


Nm 


+ 




Nt 


P3 


+ . 


(B7b) 


From Eq. (B3) the characteristic lives Np,Np 2 ,./Vp 2 , etc., can be replaced with the respective lives 
N | , /V 2 , /V^, etc., at 5 re j (orthe lives of each component that have the same probability of survival 5 re) ) 
as follows: 


In 







1 

1 

ni 

Ln 1 1 

+ 

In 

+ 


5 ref . 

5 


1 

1 


Nt, 


-tin 


+ ... 


(B8) 


where, in general, from Eq. (B3) 


and 




-tm 


N 


ref 


(B9a) 


Nt 


PI 


in 


In 


>ref 


-tin 


Nt 


, etc. 


(B9b) 
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Factoring out In (lAS re f) from Eq. (B8) gives 


N n . 


ref 


1 

m 

. i 

m 

i 

til 

1 

K 

+ 

i 

r? ■ 
1 

+ 

i 

N 3 



Mm 


or rewriting Eq. (BIO) results in 


(BIO) 


r i 

rn m 

r , i 

i 

_~N_ 

H.il 

il 

i 

4? - 
1 


(Bll) 


Substituting e as the symbol for the Weibull slope m in Eq. (Bll) gives Eq. (9) of the main text. 
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TABLE 1— TURBINE DISK MATERIAL PARAMETERS 


Material 

Titanium (Ti-6A1-4V) alloy 

Elastic modulus GPa (ksi) 

11.0(16x10') 

Poisson's ratio 

0.33 

Weibull slope 

2 

Density kg/m 3 (lb/in. ; ) 

4429 (0.16) 

Stress-life exponent 

9.2 

I 4S fatigue limit. MPa ( ksi) 

2585 (37.5) 


TABLE 2 — SUMMARY OF LOW-CYCLE-FATIGUE ENDURANCE DATA FOR 
GAS TURBINE ENGINE COMPRESSOR DISKS SUBJECTED TO BOTH 
ENGINE OPERATION AND SPIN PIT TESTS 
[Material, titanium (Ti-6AI-4V) alloy; speed, 1 1 200 rpm.] 


Test 

disk 

Operating time, cycles a 

Comments h 

Experimental life, 
cycles' 

Weibull 

slope, 

e 

Loaded 

(engine 

operation) 

Unloaded 
(spin pit 
tests) 

Total 



L 

in 

L 

so 

! Disks A (Fig. 1(a)) j 

Al 

1966 

7 100 

9 066 

Failed 

594 

7900 

219(H) 

urn 

A2 

2932 

1 I 1 00 

14 032 

Failed 









A3 

2779 

17 100 

19 879 

Failed 









A4 

3750 

21 000 

24 750 

Unfailed 









AS 

2947 

33 100 

36 047 

Failed 









Disks B (Fig. 1(b)) ! 

B1 

3437 

12 100 

15 537 

Failed 

5832 

1 3600 

19200 

5. 

5 

B2 

2616 

13 100 

15716 

Failed 









B3 

3149 

16 100 

19 249 

Failed 









B4 

2880 ! 

21 (XX) 

23 880 

Unfailed 









B5 

2755 

30 (XX) 

32 885 

Un failed 

i 









'Data from Ref. 22. 

h AU failures were bolt hole cracks exceeding 0.7938 mm (0.03125 in.) (ref. 22). 
" See Fig. 3 for Weibull plot of data 


NAS A/TM— 1 999- 1 07436 


20 




TABLK 3. — NORMALIZHD PREDICTED L u> LIVHS OF GAS TURBINE HNGINH 
COMPRESSOR DISKS AS FUNCTION LOAD 
[Material, titanium (Ti-6A1-4V) alloy; reference stress. t 4V 509.2 MPa (73.846 ksi); 
reference volume. l 7 , t . 2.264709x 10 111 m T 1 .3820 105 x 10 ' in. 1 ): reference life, 

/. |n , I cycle (normalized): disk speed, I I 200 rpm: stress-life exponent. 9.2: 
Weibull slope. 2; t 4S , fatigue limit. 258.6 MPa (37.5 ksi).j 


(a) No fatigue limit 


Parameter 

Predicted normalized L ]n life, cycles j 


■mnTi ii JITIW 

Unloaded | Loaded 

Unloaded 1 Loaded | 

Segment angle, deg | 

60 


Segment life 
Disk life 

Time at load, percent 

mm 

0.14618 

0.05968 

18 

75.118 
2 1 .7049 
85 

0.15307 

0.04419 

15 

Combined load-cycle disk life 

0.33147 

0.29123 


lb) Fatigue limit 


Segment life 
Disk life 

Time at load, percent 

8171.2 

3335.8 

82 

0.14673 

0.05990 

18 

7692.5 

2220.6 
85 

0. 15324 
0.04424 
15 

Combined load-cycle disk life 

0.33275 

0.29490 


TABLE 4.— PREDICTED L u , LIVHS OF GAS TURBINE ENGINE 
COMPRESSOR DISKS FOR ENGINE LOADS AS FUNCTION OF SPEED 
[Material, titanium (Tj- 6A1-4V) alloy; reference stress, t 4S . 509.2 MPa 

(73.846 ksi); reference volume, V rv , . 2.264709x 10 ,H m 1 < 1 .3820K)5x 1(U in. 1 ); 
stress-life exponent, 9.2; Weibull slope. 2; fatigue limit, none assumed.) 


(a) = 23833 cycles 


Speed, rpm 

Predicted life 

Disks A 

Disks B | 

Cycles 

Hi 

Cycles 

m 

3 730 


5 600 

24 153 

■ ■ 

7 459 


1 396 

2 133 

1 

1 1 200 

139 

202 

103 


14 930 

16 

23 

5 

■ 


(b) L h , = 132 373 cycles 


3 730 

21 450 

31 103 

134 150 

206 593 

7 459 

5 349 

7 754 

1 1 847 

18 245 

1 1 200 

774 

1 122 

572 

883 

14 930 

89 

128 

28 

44 


TABLE 5. — SUMMARY AND COMPARISON OF DISK L lt , LIFE PREDICTIONS WITH 
EXPERIMENTAL L ul LIVES 

[Material, titanium (Ti-6A1-4V) alloy; speed, 1 1 200 rpm.) 


Disk 

designation 




Disk L 

i life 





NAPC crack growth 

Modified crack growth 

Za ret sky probable cause 

Experimental 


method (from Ref. 22) 

method (Hq. (20)) 

method (from Table 4(b)) 

(from Table 2) 


Cycles 

Hr 

Cycles 

Hr 

Cycles 

Hr 

Cycles 

Hr 

A 


3100 


■SI 

774 

1122 

594 

861 

B 

HraS 

2400 

KS9 

■H 

572 

883 

5832 

8981 
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Speed, rpm 


16.6 9.3 



0.71 in. - J L 


(a) Disk A. 



[rjM8.6 10.5 
UI in. in. 


(b) Disk B. 

igure 1 . — Schematic of compressor disks used for 
analysis and testing. Material, titanium (Ti-6AI-4V) 
alloy. 


L 10 = 7.9x1 0 3 
L 50 = 21 .9x1 0 3 
e = 1.8 


10 20 30 50 10C 

Life, stress cycles 


(a) Disks A. 


1 2x1 0 3 


11 200 



Time, min 

gure 2. — Spin pit test cycle speed for compressor 
disks as function of time. (From ref. 22.) 


« / L 10 = 1 3.6x1 0 3 

4 “ / L so = 19.2X10 3 

/ e = 5.5 

2 I I ilili/ 1 i I i I i 1 i 

3 5 1 0 20 30 50 10C 

Life, stress cycles 
(b) Disks B. 

Figure 3. — Weibull analysis of low-cycle-fatigue 
ing of combined rim-loaded and unloaded turt 
disks. Speed, 11 200 rpm. (Data from ref. 22.) 





Normalized relative Normalized relative 

disk life, Z-io disk life, Lio 



Material Weibull slope, e 
(a) Effect of Weibull slope. 



Stress-life exponent, c 
(b) Effect of stress-life exponent. 


Figure 6. — Effect of Weibull slope and stress-life 
exponent on predicted Lio lif e of generic rotating 
disk. (From ref. 23.) 




(b) Disk B, 30° segment. 


Figure 7. — Normalized probability of failure based 
on Probable Cause analysis of rim-loaded 
compressor disks. Speed, 1 1 200 rpm. 
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Disk speed, rpm 



Figure 8. — Calculated L 0 .i lives of gas turbine engine 
compressor disks as function of speed under 
engine load. Material, titanium (Ti-6AI-4V) alloy; 
reference life, L re f, 132 373 cycles; reference stress, 
t 45 , 509.2 MPa (73.846 ksi); reference volume, V re1 , 
2.264709x1 (T 10 m 3 (1 .38201 05x1 0" 5 in 3 ); stress-life 
exponent, c, 9.2; Weibull slope, e, 2; fatigue limit, 
none assumed. 



In N 


Figure 9. — Sketch of Weibull plot where (Weibull) slope 
or tangent of line is m. Sp is probability of survival 
of 36.8 percent at which N = A/p or N/N$ = 1 . 



In N 

Figure 1 0. — Sketch of multiple Weibull plots where each numbered 
plot represents cumulative distribution of each component in system 
and system Weibull plot represents combined distribution of plots 
1 , 2, 3, etc. (All plots are assumed to have same Weibull slope m.) 
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